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13  ASS'RAC  M*» *:-a’ 

Tha  stability  characteristics  of  a  shaped-charge  jet  which  carries  an  electrical 
current  are  investigated.  The  jet  Is  assumed  to  be  perfectly  plastic,  perfectly  ! 
conducting,  and  axisymmetric .  The  equations  which  govern  the  behavior  of  the  jet  are 
solved  numerically  using  first-order  perturbation  theory,  after  the  surface  of  the  ! 
jet  has  been  subjected  to  some  arbitrary  disturbance.  The  solution  of  these  equations 
indicates  that  the  jet  is  unstable.  The  growth  rate  of  the  instability  depends  upon 
the  relative  importance  of  inertial,  electrical,  and  plastic  forces,  and  upon 
size  of  the  disturbance  wavelength  relative  to  the  jet  radius.  It  is  argued  and 
demonstrated  from  the  numerical  solutions  that  the  electromagnetic  forces  always  make 
the  jet  more  unstable  than  it  would  be  In  their  absence.  Application  of  the  theory 
suggests  that  electromagnetic  fields  may  be  capable  of  enhancing  jet  breakup  on  a 
time  scale  of  practical  interest,  but  more  work,  both  theoretical  and  experimental, 
needs  to  be  done  before  definite  conclusions  can  be  drawn.  / 
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I.  INTRODUCTION 

There  has  been  speculation,  since  Walker’s  initial  work1  on  the  problem,  that  electro¬ 
magnetic  fields  might  be  used  effectively  to  destroy  shaped-charge  jets.  Clearly  if  such  an 
effect  could  be  demonstrated  conclusively,  it  would  have  enormous  potential  for  military 
applications.  The  idea  originally  proposed  by  Walker  was  that  large  electric  currents  be 
used  to  vaporize  the  jet.  It  appears,  however,  that  electromagnetic  fields  could  be  used 
in  other  ways  as  well,  and  it  is  important  to  assess  the  effectiveness  of  various  physical 
mechanisms  contributing  to  the  destruction  of  the  jet.  In  this  reporv.  we  will  be  concerned, 
with  one  such  mechanism,  namely,  how  electromagnetic  effects  can  increase  the  growth 
rate  of  instabilities  which  have  been  shown  to  occur  even  in  the  absence  of  these  effects. 

Metallic  jets  formed  from  explosive  charges  have  been  studied  extensively.  The  basic 
theory  governing  their  formation  and  propagation  was  developed  many  years  ago  in  the 
well-known  work  of  Birkhoff,  MacDougall,  Pugh,  and  Taylor2  and,  later,  of  Pugh,  Eichel- 
berger,  and  Rostoker.3  At  least  for  a  time  after  their  formation,  the  jets  are  found  to 
stretch  as  they  propagate  with  a  reasonably  uniform  velocity  gradient  across  the  length  of 
the  jet.  Such  a  characteristic  is  highly  desirable  since  the  penetration  capability  of  the  jet 
increases  with  increasing  length.  Eventually,  however,  as  radiographs  have  demonstrated, 
the  jet  will  tend  to  neck  down  in  places  along  its  surface  and  finally  break  into  segments. 

This  behavior  suggested  that  the  idealized  stretching  motion  described  above  was 
net  stable  and  has  prompted  a  number  of  investigations  into  the  stability  characteristics 
of  jets.  Significant  among  these  were  the  early  one-dimensional  investigations  of  Chou 
and  Carleone,4  who  concluded  that  perfectly  plastic  (no  strain  hardening)  jets  should  be 
unconditionally  unstable  to  perturbations  along  their  surface.  They  identified  a  number 
of  parameters  which  affected  the  growth  rate  of  the  instability  and  compared  their  re¬ 
sults  with  computer  simulations  carried  out  by  them  and,  earlier,  by  Karpp  and  Simon.5 
These  simulations  were  undertaken  using  Wilkins’  two-dimensional,  elastic-plastic  HEMP 
code,8  Many  effects  seen  in  the  simulations  were  predicted  by  the  one-dimensional  sta¬ 
bility  analysis.  However,  the  simulations  did  suggest  the  existence  of  a  critical  or  most 
unstable  wavelength  for  the  surface  perturbation,  as  had  also  been  supposed  from  analysis 
of  the  experimental  data,  and  the  one-dimensional  stability  calculations  did  not  predict 
this  effect.  Presumably,  this  discrepancy  led  Carleone,  Chou,  and  Ciccarelli7  to  modify 
their  model  and  account  approximately  for  stress  enhancement  in  the  vicinity  of  the  neck, 
an  effect  neglected  in  their  earlier  calculations.  The  revised  version  did  predict  a  criti¬ 
cal  wavelength,  although  the  value  was  not  always  in  agreement  with  that  observed  in 
the  simulations.  Additional  one-dimensional  calculations,  along  similar  lines,  were  later 
performed  by  Walsh.8  Using  dimensioned  analysis,  he  identified  a  single  parameter  upon 
which  the  growth  rate  of  the  perturbation  depended,  and  studied  the  time  evolvement  of 
jets  with  random  initial  surface  disturbances. 

In  more  recent  work  stability  calculations  have  been  extended  to  two  dimensions  by 
Frankel  and  Weihs,9  by  Curtis,10  by  Pack,11  and  by  Romero.12  The  various  investigations 
are  all  based  on  the  assumption  that  the  jet  is  axisymmetric  and  sire  differentiated  primarily 
by  the  form  of  the  constitutive  relation  assumed.  The  most  satisfactory  assumption  seems 
to  be  that  made  by  both  Pack  and  Romero,  in  which  the  jet  was  taken  to  be  perfectly 


plastic  and  to  satisfy  the  Levy-von  Mises  stress-strain  relations.  The  work  of  Romero 
was  particularly  thorough.  He  identified  a  time-dependent  parameter  which  measured 
essentially  the  ratio  of  inertial  to  plastic  forces  in  the  jet,  and  studied  jet  stability  as  a 
function  of  both  that  parameter  as  well  as  the  structure  of  the  perturbation.  Most  of  the 
calculations  were  undertaken  via  numerical  solutions  of  the  linear  differential  equations 
which  govern  the  time  evolution  of  the  perturbation,  but  some  approximate  limiting-case 
analytic  solutions  were  also  presented. 

In  this  report  we  will  follow  basically  the  same  general  procedure  used  by  Romero. 
We  assume  a  perfectly  plastic,  shaped-charge  jet  which  initially  is  stretching  uniformly  in 
the  axial  direction.  Unlike  in  any  previous  work,  however,  we  also  assume  that  an  electrical 
current  is  conducted  along  the  surface  of  the  jet.  At  time  t  =  0  the  surface  is  subjected 
to  a  small  perturbation  so  that  its  radius  is  no  longer  independent  of  the  axial  coordinate 
z.  The  equations  of  motion  of  the  jet  are  then  solved  to  determine  how  this  perturbation 
evolves  in  time.  As  in  previous  work,  we  assume  not  only  that  the  jet  is  axisymmetric 
initially,  but  restrict  ourselves  to  perturbations  which  do  not  destroy  the  axial  symmetry. 
Our  principal  interest  lies  in  determining  how  electromagnetic  fields  associated  with  the 
current  in  the  jet  affect  the  growth  times  of  instabilities  which  occur  even  in  the  absence 
of  these  electromagnetic  effects. 

A  number  of  other  assumptions  are  made  in  this  first-effort  calculation  in  order  to 
make  the  analysis  as  tractable  as  possible.  First,  we  assume  as  in  Refs.  9  and  12  that  the 
jet  is  infinitely  long.  In  practice,  this  assumption  means  only  that  the  dynamics  of  the  tip 
and  tail  of  the  jet  can  be  ignored,  and  implies  that  the  analysis  holds  only  at  points  well 
removed  from  the  ends.  This  assumption  seems  reasonable  because  the  length  of  the  jet 
is  large  compared  *o  its  radius  and,  as  will  be  seen,  disturbances  which  grow  appreciably 
have  initial  wavelengths  which  are  smaller  than  or  comparable  to  the  initial  radius  the  jet. 
Second,  we  have  assumed  that  the  jet  is  a  perfect  conductor.  This  assumption  is  strictly 
valid  only  if  the  magnetic  Reynolds  number13  associated  with  the  jet  is  sufficiently  large, 
a  condition  that  is  questionable  for  the  case  at  hand.  It  is  common,  however,  in  stability 
calculations  to  assume  infinite  conductivity  anyway,  perhaps  in  the  belief  that  the  basic 
conclusions  reached  will  not  be  significantly  different  from  those  obtained  in  the  finite- 
conductivity  case.  In  assuming  infinite  conductivity,  we  are  in  effect  replacing  J  x  B  body 
forces  by  surface  forces  which  result  from  current  being  carried  only  on  the  surface  of  the 
jet.  We  are  also  neglecting  increases  in  the  internal  energy  of  the  jet  result:ng  from  resistive 
heating.  Consequently,  we  can  expect  the  analysis  to  hold  only  for  currents  sufficiently  low 
and  times  sufficiently  short  that  large-scale  changes  in  the  energy,  such  as  would  produce 
melting,  do  not  occur.  These  points  will  be  further  discussed  in  Secs.  V  and  VI.  Finally, 
we  will  neglect  dissipation  by  plastic  forces  and  assume  that  the  jet  is  isothermal.  Thus, 
we  will  not  consider  any  energy-conservation  equation  in  our  analysis  and  this  omission 
also  facilitates  the  calculations. 

The  organization  of  the  report  is  a s  follows.  In  Sec.II,  we  describe  the  model  and 
indicate  the  basic  equations  assumed  to  govern  the  behavior  of  the  jet.  In  Sec.  Ill,  the 
lowest-order  solution  to  these  equations  is  worked  out.  In  Secs.  IV  and  V,  the  first-order 
equations  are  derived  and  solved,  and  the  stability  characteristics  of  the  jet  discussed. 
Finally,  Sec.  VI  contains  our  conclusions  and  plans  for  future  work. 
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These  equation*  have  the  following  meaning*.  Equation*  (1)  (3)  expreaa  con*erv»tion 
of  momentum  and  ma*«.  The  jet  hu  been  aasumed  to  be  mcompre**jble  no  that  mm 
conservation  take*  the  simple  form  represented  by  (31.  The  last  term  on  the  ri*ht-hand 

of  )>oth  EiqaU)  and  <2)  iepre*ent«  the  Lorentx  force.  The  jet  1*  uiumed  to  be  electrically 

neutral  so  that  electrostatic  force*  need  not  be  included  in  theae  equation*  and,  in  fact, 
any  other  external  force*  are  taken  to  be  ne^lijible.  Equation*  (4)-(7)  are  the  *tandard 
equation*  which  govern  the  motion  of  a  perfectly  plaatic  axi»ytnxnetnc  material  in  cylin- 
inrai  coordinate*,  with  E^q.  (4)  (0)  representing  the  Levy-von  Mi*e»  *tre«*- strain  relation* 
and  Eq.(7)  the  von  Mine*  yield  condition.  The  parameter  v  is  a  time- dependent  factor 
which  must  (>e  determined  in  the  Hnalysi*.  Equation*  (8)-(  10)  are  Maxwell’s  equation*. 
H  e  displacement  current  ha*  twen  neglected  in  £q.  (10);  in  *eneral,  mcludin*  this  term 
produce*  only  terms  of  the  order  of  c/c,  where  i>  is  a  typical  velocity  in  the  jet  and  c  i*  the 
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light  speed.  Consequently,  the  effect  of  this  current  is  clearly  negligible  for  the  problem  at 
hand.  Equations  (11)  and  (12)  represent  Ohm’s  law  for  the  two  nonvanishing  components 
of  the  current  density. 

In  the  vacuum  stirrounding  the  jet,  i.e.,  for  r  >  r»,  Eqs.  (3)-(7)  are  meaningless,  but 
the  remaining  equations  hold  with  p  =*  0,  J  »  0,  and  <raa  *  0.  These  equations  therefore 
imply  certain  boundary  conditions  that  must  exist  at  the  interface  between  the  jet  and  the 
vacuum.  The  appropriate  conditions  can  be  obtained  by  integrating  each  of  the  relevant 
equations  across  the  boundary,  of  assumed  thickness  6,  in  a  direction  normal  to  the  surface 
and  then  taking  the  limit  as  6  — *  0.  Let  n  be  the  unit  normal  to  the  separating  surface 
pointing  in  the  positive  radial  direction  and  let  dt  be  an  element  of  arc  length  normal 
to  the  surface. ,  Following  the  procedure  discussed  elsewhere,11  we  now  multiply  each  of 
the  pertinent  equations  through  by  dt  and  integrate  across  the  boundary.  We  allow  for  a 
surface  current  j*  so  that  the  integral  of  J  produces  j*  as  6  — *  0.  Finally,  we  observe  that 
the  integral  of  each  partial  time-derivative,  ^j-,  produces,  except  for  terms  which  vanish 
with  6,  the  result  - v  ■  n  <  F  >.  The  brackets  here  denote  the  change  in  F  as  the  boundary 
is  crossed. 

The  following  results  are  obtained  when  the  integrations  above  are  carried  out.  After 
we  make  use  of  Eqs. (8)  and  (9),  the  integrals  of  Eqs.  (1)  and  (2)  yield 
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Equations  (11)  and  (12)  yield  ;*-=0  unlr-!*  rr  is  infinite  in  which  case  they  yield  no  infor¬ 
mation.  We  finally  observe  that  since  the  boundary  moves  with  the  jet,  we  must  have 
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and  that  the  unit  normal  n  is  given  by 
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This  system  of  twelve  equations  and  appropriate  boundary  conditions  are  sufficient 
in  principle  to  completely  specify  the  motion  of  the  jet  or,  alternatively,  to  determine  aa  a 
function  of  position  and  time  the  twelve  unknowns:  u,  w,  9r,  9 ,,  9t,  9,,,  v,  E,,  E,,  Jrt  J », 
and  B.  The  equations  are  obviously  difficult  to  solve  in  general,  and  we  will  be  interested  in 
simplifying  them  by  undertaking  a  linear  perturbation  solution.  Through  such  a  solution  it 
is  possible  to  study  the  stability  of  the  motion  of  the  jet  to  various  type*  of  perturbations  in 
a  comparatively  simple  way.  Furthermore,  as  will  be  seen,  it  is  sometimes  possible  to  solve 
only  a  subset  of  the  available  equations,  the  remaining  ones  being  necessary  to  determine 
properties  of  the  jet  which  do  not  affect  the  stability  analysis. 


III.  ZERO-ORDER  SOLUTION:  IDEALIZED  MOTION  OF 

THE  JET 


Experimental  studies  in  the  motion  of  shaped-charge  jets  have  indicated  that  at  early 
times  the  jets  stretch  rather  uniformly  as  they  propagate,  with  the  jet  radius  being  in¬ 
dependent  of  z  and  the  velocity  w  increasing  linearly  across  the  length  of  the  jet.  At 
later  times  there  appear  along  the  surface  of  the  jet  perturbations  which  tend  to  disrupt 
this  idealized  motion.  The  purpose  of  stability  analysis  is  to  determine  whether  these 
perturbations  grow  in  time  and,  if  so,  how  rapidly. 

In  order  to  study  jet  stability  in  the  presence  of  electroma<pietic  fields,  we  will  assume 
the  same  basic  type  of  idealized  motion  at  early  times.  Thus  we  take  r\  independent  of 
z  and  w  *  A(t)z,  where  A(f)  is  to  be  determined.  As  indicated  previously,  it  will  also 
be  assumed  in  this  first  analysis  that  the  jet  is  perfectly  conducting  (<7  infinite)  and  that 
the  current  I  is  carried  along  the  outer  surface,  i.  e.,  at  r  w  r».  We  take  the  initial  strain 
rate  to  be  3  (assumed  given),  the  initial  radius  to  be  r o,  and  denote  by  subscript  zero  the 
lowest-order  solution  to  the  governing  equations,  i.  e.  the  solution  which  represents  the 
idealized  motion.  The  superscript  V  is  used  hereafter  to  denote  the  value  of  some  quantity 
in  the  vacuum,  whereas  the  absence  of  a  superscript  denotes  the  value  in  the  jet.  The 
origin  of  the  coordinate  system  can  be  chosen  arbitrarily  and  we  take  one  which  moves 
with  the  center  of  moss  cf  the  jet.  Thus,  w  in  the  governing  equations  means  the  velocity 
relative  to  the  center  of  mass  and  w  =  0  at  x  =»  0. 

According  to  th?  assumptions  above  and  from  Ampere's  law,  it  follows  that  in  the  jet 
Jo  =*  0  and  Bo  =*  0,  while  in  the  vacuum 


b: 


m 


Finally,  since  n,0  —  1  and  n,  0  =  0,  we  conclude  from  the  boundary  con.L  -Lons  in  (15)  and 
1 16)  that 
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Table  1.  Solution  of  Governing  Equations  in  Zero- Order. 


Jet 

Vacuum 

Surface 

U>0  =  102 

pv _  iti 

B0  ~  iTr 

nT,o  =  1 

u0  =  ~10r/2 

n*,o  =  0  1 

n.o  =  7 1/2ro 

<?r.o  -  3£iiJ~( r 2  -  7ro)  -  gife 

Jrrcrr1^3 

«...  =  -  r-8)  - 

<7„,0  =  0 

v0  =  $  ..  . 

Jr,  0  —  0 

J,.o  —  0 

O 
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a? 

and 


(21) 


r,0 


=  0. 


(22) 


The  remaining  fluid- mechanical  properties  of  the  jet  in  zero  order  follow  by  solving  the 
equations  of  the  preceding  section  under  the  assumptions  discussed  above.  The  analysis 
is  similar  to  that  which  has  been  undertaken  elsewhere11,1*  and  will  not  be  repeated  in 
detail  here.  The  only  significant  difference  with  the  nonelectrical  case  is  that  the  stress 
components  must  satisfy  the  boundary  conditions  represented  by  <T3)  and  (14),  whereas 
for  no  magnetic  fields  the  right-hand  sides  of  these  equations  reduce  to  zero.  Results  of  all 
zero-order  quantities,  except  for  the  electric  fields  which  are  not  needed  in  the  subsequent 
analysis,  are  presented  for  reference  in  Table  I.  They  can  be  seen  by  direct  substitution  to 
satisfy  the  equations  and  boundary  conditions  in  Sec.  II.  In  this  table  the  parameter  7  is 
given  by 


1 

1  +  (Sf 


(23) 


r 


IV.  FIRST-ORDER  EQUATIONS 


We  now  wish  to  consider  the  form  taken  by  the  equations  of  Sec.II  in  first-order 
perturbation  theory.  We  therefore  assume  that  each  physical  quantity  can  be  represented 
by  the  zero-order  solution  plus  a  small  correction,  and  linearize  the  governing  equations 
about  those  corrections.  Therefore,  for  each  quantity  F  we  assume 


F  zz  F0  +  Fi(r,z,t),  (24) 

substitute  in  the  equations  of  Sec.II,  and  retain  only  the  terms  which  are  of  order  F\  or 
lower.  The  solution  of  the  first-order  equations  should  then  predict  the  stability  of  the  jet 
to  the  idealized  motion  described  in  the  preceding  section. 

Considering  first  just  the  fluid-mechanical  equations  in  the  jet,  we  find  after  making 
the  substitution  above  in  Eqs.(l)-(7),  that  the  resulting  equations  can  be  partly  uncoupled 
to  produce 


du\  du0  9u  i  d<r,\  1  d*u\ 

“ST  +  aT  +  aT  +  “  HT  +  iZ IF 


(25) 


dw i  dw\  du’o  dw\  do,  i 

P"zr~  4-  puo-^r-  +  pWi  4-  pti'o—  -  —  4- 

Oi  dr  O :■  dz  dz 


1  .d^Wi  1  dw  j 
6i/0  dr*  +  r  dr 


dz *  J 


(26) 


and 


du\  Ui  dw\ 
dr  ^  r  ^  dz 


(27) 


Similarly,  the  electromagnetic  equations,  (8)-(12),  can  be  uncoupled  to  yield  a  single  result 
which  describes  the  convection  of  the  induction  field  in  the  jet,  namely, 


dB , 
dt 


du0 

dr 


B\  4-  up 


dB\ 

dr 


4-  wp 


dB  i 
dz 


du'o 

4-—  B,  =  0. 
dz 


(28) 


The  results  from  Table  I  are  to  be  used  to  obtain  the  zero-order  functions. 

In  the  vacuum  it  seems  that  there  can  be  no  perturbation  to  the  existing  induction 
field.  This  conclusion  follows  most  directly  from  Ampere’s  law  and  the  assumed  symmetry 
of  the  problem.  That  is,  since  there  is  only  an  azimuthal  component  of  B  and  this  com¬ 
ponent  is  independent  of  angle,  the  induction  field  is  given  to  all  orders  by  Eq,  (20)  and 
we  have 


B\  =  0. 


(29) 
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As  indicated  before,  we  will  find  no  need  in  this  problem  to  calculate  the  first-order  electric 
fields  in  the  vacuum,  provided  we  restrict  ourselves  to  studying  the  stability  properties  of 
the  jet. 

Finally,  we  consider  the  boundary  conditions  in  first  order.  Evidently,  we  have  from 
Eq.(19) 


«i  =  — g^a,.  (30 ) 

Other  first-order  quantities  evaluated  at  the  boundary  are,  in  accordance  with  standard 
procedure,  represented  by  their  Taylor  expansion  about  r  =  t\0-  Consequently,  for  any  F 
we  have  through  first  order 


F{rb)  =  F0(r6,0)  +  Fi(r*,0)  +  (31) 

We  now  use  Eq.(31)  and  the  results  in  Table  I  to  evaluate  the  boundary  condition  repre¬ 
sented  by  Eq.(13).  We  find 

_  /  )  _  3p/?275/2r0v  _  y  ,/^ix  _  jjiKo)} 

*'o)  l4ir273/2r3  ~  4  *  M  37/?f(  dr  )r‘°  r6,0  (32) 

In  obtaining  this  result  we  have  used  the  expression 

A 

-aT  -  ~  =  6l/O(0V,l  -  <7x,i),  (33) 

Or  r 

easily  derived  from  the  linearized  versions  of  (4), (5),  and  (7),  to  express  <7r,i  in  terms  of 
<7, ,i  and  tij.  Similarly,  from  Eq.  (14)  and  the  linearized  version  of  (6),  we  have 

37/3l  ar  )'M  +  ^  ar  ^  dz  '  ^ 

Finally,  we  have  from  (18)  in  first  order 


arM  ,  a  drhil  ,  10 

~df  +  7*z"aT  =  u'(ri0>  ~  ~2rh'1'  (35) 

The  remaining  conditions,  those  implied  in  (15)-(17),  need  not  be  considered  in  the  sub¬ 
sequent  analysis. 

We  now  observe  that  Eqs.  (25), (26),  and  (27)  along  with  the  boundary  conditions 
represented  by  (32),  (34),  and  (35)  are  sufficient  to  determine  the  quantities  tij,  w\,  <rXil, 
and  r*  !  for  all  time,  given  some  set  of  initial  conditions.  As  will  be  seen  these  results  are 
sufficient  to  study  the  stability  of  the  jet. 
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It  is  convenient  for  further  analysis  to  follow  the  basic  procedure  of  Romero12  and 
perform  a  transformation  which  removes  the  explicit  z  dependence  from  the  governing 
equations.  Consequently,  we  define  new  coordinates 


(36) 


and  substitute  into  Eqs.  (25), (26), (27), (32), (34), and  (35).  Since  the  results  now  contain 
no  explicit  dependence  on  £,  the  solution  of  the  equations  can  be  taken  to  be  of  the  form 


t)  -  0roui(ri,  r)e,Ki 

u>i(*7,s,r)  =  0rowt(T],T)e'Ki 

<r.Av,  Z,T)  =  Y&tAv*  r)e'K< 

n,i(Z,T)  =  r0r»,i(r)e,x<.  (37) 

Here  A"  is  a  dimensionless  wave  number,  related  to  the  real  wave  number  k  through  the 
relation 

K  =  r0k  =  (38) 

Ao 

where  A0  is  the  initial  wavelength  of  the  perturbation.  We  have  also  normalized  the 
velocities  by  0ro,  the  stress  by  Y,  the  jet  radius  by  r0,  and  defined  a  dimensionless  time  r 
given  by 


r  =  Bt. 


(39) 


Using  these  results  in  Eqs.(25),  (26),  (27),  (32),  (34),  and  (35),  we  find  that  the  dimen¬ 
sionless  first-order  variables  must  satisfy  the  following  set  of  equations 


du)\ 

dr 


4-  7ti>i  =  iAK~r&,' 


d3u’\ 
dr B 


+  i^L+1]A-w,|, 


+  Hi  +  =  0, 

dr,  rj 


(40) 

(41) 


and 


dPvt<\  1  <9<7x, j  ■?,, -2.  iK  ,dBw\  1  dw\  , 

- +  — - - rAVi  =  ^ - 1  1 


dr,2  r,  dr, 
with  boundary  conditions 


3  1  dr,2  +  q'dn 


+  73(2A'2-9/.4)u)1j,  (42) 
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(43) 


V.  SOLUTION  OF  FIRST-ORDER  EQUATIONS  AND 
STABILITY  CHARACTERISTICS  OF  THE  JET 

We  now  proceed  to  solve  Eqs.  (40)-(45)  of  the  preceding  section  and  thereby  to  study 
the  stability  characteristics  of  the  jet  as  a  function  of  the  dimensionless  parameters  A,  H, 
and  K.  Physically,  A  and  H  measure,  respectively,  the  importance  of  plastic  and  electrical 
forces  relative  to  inertial  forces,  and  K  is  a  measure  of  the  size  of  the  initial  wavelength  of 
the  perturbation  relative  to  the  initial  radius  of  the  jet.  It  is  noteworthy  that  the  parameter 
H,  and  thus  all  electrical  effects,  occurs  only  in  the  boundary  condition  represented  by 
Eq.(43).  This  rather  special  dependence  is  a  consequence  of  our  assumptions  of  infinite 
conductivity  and  axial  symmetry.  As  will  be  seen  in  the  following  analysis,  however,  this 
one  term  can  have  a  marked  effect. 

Before  the  equations  can  be  solved  we  must  first  specify  a  set  of  initial  conditions. 
Apparently,  one  need  only  specify  Wi(tj,t  =  0)  since  the  remaining  conditions,  namely, 
those  on  fij,  h.i,  and  cr2,i,  then  follow  from  Eqs.(41),  (44),  and  (43),  respectively.  In  the 
remainder  of  the  analysis  we  restrict  ourselves  to  initial  values  of  tZ>i  which  are  constant 
and  denoted  by  iC  where  C  is  real.  Solving  the  equations  noted  above  then  yields 
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CKr> 

0,(77,  r  =  0)  =  — — , 


(48) 


r\,(r  =  0)  =  (49) 

and 


-  0)  =  l~(H  -  1)  +  (2 K1  -  9 -  (2JC1  -  9/ A)^  (50) 

where  Iq  is  a  modified  Bessel  function.  The  constant  C  can  be  chosen  arbitrarily  and  we 
take  C  =  6/A'  so  that  f*tl(r  =  0)=1. 

To  solve  Eqs.  (40)-(45),  we  have  represented  both  time  and  space  derivatives  by  finite 
differences  and  numerically  solved  the  resulting  set  of  linear  algebraic  equations  by  matrix 
inversion.  The  procedure  is  quite  straightforward  and  need  not  be  discussed  here  in  detail. 
Typically,  a  time  step  A r  of  0.005  and  a  grid  spacing  which  depended  on  the  wavelength  of 
the  initial  perturbation  was  used.  In  general,  small  wavelengths  require  more  grid  points 
to  achieve  spatial  resolution  and  we  have,  in  some  cases,  used  as  many  as  400.  This  point 
has  been  further  discussed  by  Romero.12 

In  stability  calculations  in  which  there  is  explicit  time  dependence,  it  is  customary  to 
consider  how  the  relative  amplitude  u;  given  by 


7"6.0 


(51) 


varies  with  time.  If  this  quantity  grows  monotonically  the  jet  is  assumed  to  be  unstable, 
whereas  if  it  remains  constant  or  approaches  some  asymptotic  value,  the  jet  is  stable. 
We  therefore  wish  to  consider  this  amplitude,  obtained  by  solving  the  equations  discussed 
previously,  for  several  cases  of  interest. 

Shown  in  Fig. 2  is  as  a  function  of  r  for  different  values  of  the  dimensionless  param¬ 
eter  K.  The  parameters  A  and  H  were  taken  to  be  1.5  and  zero,  respectively.  The  value 
of  A  is  reasonable  since  this  quantity  is  of  order  unity  for  typical  shsoed-charge  jets  and 
H  —  0  corresponds,  of  course,  to  no  electrical  current  present.  For  typical  strain  rates, 
each  urut  of  dimensionless  time  corresponds  in  real  time  to  a  few  tens  of  microseconds. 
As  is  evident  for  each  value  of  K  considered,  u  continues  to  grow  for  the  time  plotted. 
Consequently,  the  jet  is  unstable  regardless  of  the  initial  wavelength  of  the  perturbation. 
The  value  K  —  ir,  which  corresponds  to  a  wavelength  equal  to  the  initial  diameter  of  the 
jet,  leads  to  the  maximum  initial  growth  rate  for  the  various  cases  considered.  The  two 
shortest  wavelengths,  namely,  those  corresponding  to  K  =  2x  and  K  =  4t  are  nearly 
stable  initially.  However,  their  amplitudes  begin  to  grew  once  the  jet  has  stretched  enough 
that  the  wavelength  becomes  comparable  in  size  to  the  jet  diameter.  At  that  time  the 
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amplitude  grows  very  rapidly,  eventually  becoming  larger  them  the  amplitude  of  the  longer- 
wavelength  perturbations  which  had  high  initial  growth  rates.  For  the  two  long- wavelength 
modes  considered,  corresponding  to  K  given  by  t/2  and  tt/8,  the  growth  rate  decreases 
with  increasing  wavelength.  These  results  are  in  basic  agreement  with  the  conclusions 
reached  by  Romero.12 
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Figure  2.  Relative  amplitudes  versus  time  for  different  values  of  K  and  with  no  electrical 
current  present.  The  parameters  A  and  H  were  given  by  1.5  and  zero,  respectively. 

The  corresponding  cases  for  H  ^  0,  i.e.,  for  an  electrical  current  present,  are  shown 
in  Fig.  3.  For  currents  of  the  order  of  100  kA  or  so,  the  parameter  H  is  also  of  the  order 
of  unity  for  typical  jets.  For  the  cases  plotted  in  the  figure,  H  was  actually  taken  to  be 
somewhat  smaller,  namely,  0.27.  Apparently,  the  current  produces  no  change  in  the  basic 
shapes  of  the  curves  or  in  the  qualitative  conclusions  reached  for  the  no-current  case.  It  is 
evident,  however,  that  the  current  does  produce  an  increase  in  the  relative  amplitude  at 
any  given  time  for  all  the  wavelengths  considered.  This  effect  can  be  seen  most  easily  in 
Fig.4  in  which  the  ratio  w/fjio/u«=o,  with  appropriate  values  of  u  taken  from  Figs.  2  and  3, 
is  plotted  for  each  value  of  K.  Clearly  this  ratio  corresponds  to  the  ratio  of  the  perturbed 
radius  amplitude,  f6ii,  in  the  current  and  no-current  cases.  Evidently,  the  overall  effect  of 
the  electrical  current  is  to  enhance  the  instability. 

Our  qualitative  understanding  of  this  instability  and  of  the  behavior  exhibited  in  Figs. 
2  and  3  can  be  described  approximately  as  follows.  When  the  jet  begins  to  neck  down  in 
some  localized  region,  there  is  a  tendency  for  the  material  on  the  right-hand  side  of  the 
neck  minimum  to  be  accelerated  toward  regions  of  increasing  cross  sectional  area,  and  to 
be  decelerated  on  the  opposite  side.  If  there  were  no  change  in  the  stress,  then,  the  mass 
in  the  neck  would  be  decreased  still  more,  resulting  in  unconditional  instability.  In  reality, 
however,  the  longitudinal  stress  in  the  vicinity  of  the  neck  is,  on  the  average,  increased  by  an 


13 


Figure  3.  Relative  amplitude  u>  versus  time  for  different  values  of  K  and  with  an  electrical 
current  present.  The  parameters  A  and  H  were  given  by  1.5  and  0.27,  respectively. 


Figure  4.  Ratio  of  perturbed  radius  amplitudes  for  tue  current  and  no-current  cases.  The 
parameter  A  was  given  by  1.5 
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amount  which  increases  with  increasing  curvature  of  the  surface  or  decreasing  wavelength  of 
the  perturbation.  Since  the  material  prefers  to  flow  from  regions  of  low  stress  to  regions  of 
high  stress,  this  enhancement  tends  to  cause  flow  back  into  neck  and  produce  a  stabilizing 
effect.  For  sufficiently  short  wavelengths,  near  stability  can  result.  The  wavelength  does, 
however,  increase  with  time  because  the  material  is  stretching  and  perturbations  which 
are  almost  stable  initially  eventually  become  very  unstable.  For  extremely  long- wavelength 
disturbances, the  perturbation  is  hardly  felt  in  local  regions  of  the  jet.  Consequently,  the 
growth  rate  of  the  instability  decreases  with  increasing  wavelength,  with  the  jet  eventually 
becoming  neutral  to  all  perturbations  in  the  infinite  wavelength  limit.  When  a  current  is 
added,  the  current  must  flow  down  the  left-  hand  side  of  the  neck  and  back  up  the  right- 
hand  side.  The  axial  component  of  the  Lorentz  force  associated  with  this  current  points 
in  the  positive  z  direction  on  the  the  right  and  in  the  negative  z  direction  on  the  left.  The 
effect  of  the  current  therefore  is  to  accelerate  material  out  the  neck,  always  increasing  the 
perturbation  and  enhancing  the  instability. 

In  practical  applications  of  electromagnetic  forces  to  shaped-charge  jets,  it  is  likely 
that  the  current  will  be  applied  at  some  distance  from  where  the  jet  is  initially  formed. 
Obviously  the  point  of  these  calculations  has  not  been  to  simulate  realistically  the  details 
of  such  an  occurrence  but,  rather,  to  examine  how  the  current  affects  the  general  stability 
characteristics  of  the  jet.  Some  estimate  of  the  response  of  the  jet  subsequent  to  current 
initiation  can  be  obtained  from  the  foregoing  theory,  however,  by  taking  for  the  initial 
strain  rate  0  and  initial  radius  r0  values  appropriate  for  the  jet  when  the  current  is  initiated. 
These  values  are  somewhat  smaller  than  the  real  initial  values.  It  is  then  implicitly  assumed 
that  prior  to  current  initiation  the  surface  of  the  jet  is  negligibly  disturbed. 

To  undertake  such  a  calculation  we  considered  a  copper  jet  similar  to  that  discussed 
in  Ref.  11,  and  assumed  that  the  current  was  applied  when  the  radius  of  the  jet  was 
2  mm.  A  simple  calculation  then  revealed  that  the  strain  rate  at  that  time  was  about 
2.5  x  104s_I.  Other  parameters  taken  in  the  calculation  were  p  =  8.9  x  103  kg/m3  and 
Y  —  100  MPa,  and  a  current  /  of  150  kA  was  assumed.  We  then  found  from  Eqs.  (46)  and 
(47)  that  A  ~  4  5  and  H  ~  10.7.  Results  of  the  numerical  solution  of  Equations  (40)-(45) 
for  these  parameters  are  shown  in  Fig. 5  where  u>  is  plotted  versus  the  real  time  t  for  both 
the  current  and  no  current  (H=0)  cases.  The  dimensionless  wavelength  K  was  taken  to 
be  3jt/4,  a  value  reasonably  close  to  the  most  unstable  one.  The  electromagnetic  forces 
obviously  produce  quite  a  remarkable  increase  (more  than  a  factor  of  two)  in  disturbance 
amplitude  in  the  fairly  short  time  of  about  25 fis.  It  seems  reasonable  to  expect  therefore 
that  these  forces  could  contribute  significantly  to  the  breakup  time  of  the  jet.  We  should 
also  point  out  the  the  assumed  current  is  fairly  modest  and  is  not  expected  to  produce 
large-scale  heating  of  the  jet  during  this  time  period. 


VI.  SUMMARY,  CONCLUSIONS,  AND  FUTURE  WORK 

We  have  developed  a  model  and  undertaken  calculations  to  examine  the  effects  of 
electromagnetic  forces  on  the  stability  of  a  perfectly  conducting  shaped-charge  jet.  Basic 
conclusions  reached  as  a  result  of  this  study  can  be  summarized  as  follows:  (l).In  general, 


▼ 

i 


15 


t  (ms) 


Figure  5.  Relative  amplitude  u>-  versus  time  for  current  and  no-current  cases.  The  pa¬ 
rameters  .4  and  I\  were  given  by  4.5  and  3rr/4,  respectively. 

the  jets  become  more  unstable  as  either  of  the  dimensionless  parameters  .4  or  H  increases. 
The  first  parameter  is  a  measure  of  the  magnitude  of  plastic  forces  relative  to  inertial 
forces,  and  the  second  the  magnitude  of  electrical  forces  relative  to  inertial  forces.  (2). 
The  stability  of  the  jet  also  depends  on  how  the  wavelength  of  the  perturbation  compares 
in  size  relative  to  the  radius  of  the  jet.  Both  the  very  short  and  the  very  long  wavelengths 
are  more  stable  than  those  of  moderate  magnitude,  and  there  is  a  "critical”  wavelength 
which  produces  maximum  instability.  (3).  The  stability  properties  of  the  jet  vary  with  time 
because  the  jet  stretches  as  it  propagates.  Thus,  for  example,  as  the  jet  radius  decreases 
with  increasing  time,  electrical  effects  become  more  and  more  important  because  of  the 
increased  induction  field  at  the  surface  of  the  jet.  (4).  Many  factors  make  the  suitability  of 
the  calculations  tenuous  for  computing  in  detail  how  electromagnetic  fields  will  affect  jet 
breakup  time.  However,  what  are  thought  to  be  reasonable  estimates  for  the  parameters 
involved  lead  to  predictions  that  these  forces  should  indeed  produce  a  significant  effect  on 
a  time  scale  of  practical  interest. 


In  the  future  it  would  be  of  interest  to  the  extend  the  above  calculations  to  cases 
which  do  not  preserve  the  axial  symmetry,  and  to  examine  the  effects  of  finite  conductivity 
on  the  growth  rates  of  the  instabilities.  The  first  of  these  problems  is  currently  under 
investigation  in  house  and  the  second  is  slated  for  future  study.  Obviously,  it  would  also 
be  beneficial  to  have  some  experimental  data  on  jet  disruption  to  compare  with  the  theory. 
Some  experimental  work16  has  already  been  carried  out  and  it  is  hoped  that  more  will  be 
forthcoming  in  the  near  future. 
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